Dephasing enhanced transport in non-equilibrium strongly-correlated quantum systems 
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A key insight from recent studies is that noise, such as dephasing, can improve the efficiency of quantum 
transport by suppressing coherent single-particle interference effects. However, it is not yet clear whether de- 
phasing can enhance transport in an interacting many-body system. Here we address this question by analysing 
the transport properties of a boundary driven spinless fermion chain with nearest-neighbour interactions sub- 
ject to bulk dephasing. The many-body non-equilibrium stationary state is determined using large scale matrix 
product simulations of the corresponding quantum master equation. We find dephasing enhanced transport only 
in the strongly interacting regime, where it is shown to induce incoherent transitions bridging the gap between 
bound dark-states and bands of mobile eigenstates. The generic nature of the effect is illustrated and shown not 
to depend on the integrability of the model considered. As a result dephasing enhanced transport is expected to 
persist in more realistic driven systems of strongly-correlated particles. 



Introduction — Recently the effects of noise on the effi- 
ciency of quantum transport phenomena have been scrutinised 
intensely by the scientific community. This has been moti- 
vated in part by a series of ground-breaking nonlinear spec- 
troscopic experiments on light-harvesting complexes demon- 
strating surprisingly long-lived quantum coherence during ex- 
citon transport, even in a warm and wet environment [1-3]. 
Yet transport for purely coherent exciton dynamics in such 
protein pigment networks is highly suppressed due to destruc- 
tive interference between different propagation pathways. In- 
stead studies revealed that the remarkably high transport effi- 
ciency observed (above 95%) in fact emerges in combination 
with local noise, such as dephasing, which disrupts this inter- 
ference opening up previously inhibited pathways for trans- 
mission [4—7]. Transport properties in such open systems thus 
not only defy the traditional understanding of when quantum 
effects should play a significant role, but also challenge the 
notion that couplings to the environment unconditionally de- 
grade quantum phenomena. 

Different effects of dephasing have been studied in net- 
works populated by single particles, in scenarios such as trans- 
port through quantum optical systems [8], information trans- 
mission [9] and quantum information processing [10]. How- 
ever, the interplay between noise and strong correlations in- 
duced by interactions in a many-body setting is not yet fully 
understood. In particular a recent experiment [11] with cold- 
atoms exploring these effects during time-of-flight expansion 
raises an important question as to whether and how dephasing 
can enhance transport in an interacting system. 

Here we answer this question in the affirmative by con- 
sidering a concrete example composed of spinless fermions 
with nearest-neighbour interactions hopping through a tight- 
binding chain, as depicted in Fig. 6. This model makes an 
ideal testbed for several reasons. First, it is equivalent to the 
well studiedXXZ spin-1 /2 chain [12-14], representing one of 
the simplest models of strongly-correlated electron systems, 
and its physics is still not fully understood despite being Bethe 
ansatz solvable [15]. Second, the transport properties of such 
low-dimensional interacting quantum systems remains an im- 
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FIG. 1: Spinless fermions hop with amplitude T across a chain, sub- 
ject to a nearest-neighbour density-density interaction of strength A, 
local dephasing at a rate y, and boundary driving that injects/ejects 
fermions at a rate proportional to F and driving bias /. 



portant open problem, exhibiting anomalous features such as 
ballistic spin transport and unusually high thermal conduc- 
tivity, reported experimentally in so-called spin-chain mate- 
rials [16-19]. In addition to solid state systems like chains 
of coupled quantum dots [20] or molecular wires embedded 
between electrodes [21], understanding this model is directly 
relevant to ion-trap [24], coupled-cavity array [25], and cold- 
atom [11, 22, 23] quantum systems. Of particular relevance 
are recent seminal experiments that revealed contact and bulk 
resistivity of cold fermionic atoms flowing within a narrow 
mesoscopic channel between a pair of reservoirs [26, 27]. 

Very similar to this cold-atom experiment we consider 
a chain attached to two unequal Markovian reservoirs 
at its boundaries providing continuous incoherent driving, 
along with local dephasing noise along its extension (see 
Fig. 6) [28-31]. The transport properties of the chain are 
then found by computing the current-driving characteris- 
tics of its non-equilibrium stationary state (NESS). Being a 
one-dimensional and homogeneous chain, previously studied 
single-particle interference effects originating from geometry 
or disorder are absent [5, 6, 32]. Surprisingly, we show that 
dephasing-enhanced transport nonetheless emerges so long as 
the interactions are strong enough. In this regime the NESS 
forms a cooperative many-body quantum state with a long- 
ranged domain of particles pinned to one boundary strongly 
suppressing the current, analogous to a Coulomb or Pauli 
blockade insulator [29]. Even a small dephasing rate in- 
duces incoherent transitions between this bound state, prefer- 
entially occupied by the driving, and bands of mobile scatter- 



2 



ing states. This effect significantly enhances transport, most 
especially for regimes where the current is otherwise vanish- 
ingly small. We isolate the essential conditions for enhance- 
ment and demonstrate its generality beyond the integrable sys- 
tem considered. Furthermore, the mechanism described helps 
explain recent experimental results with cold-atoms where the 
expansion of a strongly-interacting atomic gas was found to be 
slow in the absence of noise and substantially increased once 
noise was added [11]. 

Model — We study the interacting spinless fermion chain 
described by the Hamiltonian (taking /i = 1 throughout) 
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where cj , Cj are standard fermionic creation/annihilation op- 
erators for site j, nj = CjCj is the associated number opera- 
tor, and is the number of sites. In addition to the tight- 
binding hopping terms, with amplitude T, this Hamiltonian 
has a nearest-neighbour density-density interaction 1 
which, depending on the sign of its strength A, either penal- 
izes (A > 0) or promotes (A < 0) the adjacency of fermions. 
We take X = 1 to set the energy scale to which all other quan- 
tities are measured. The dynamics of the system is described 
by a Lindblad quantum master equation [33] 
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where p is the density matrix of the chain, and L is the dis- 
sipator describing the coupling to the Markovian reservoirs. 
In Lindblad form the dissipator can be expressed as L{p) = 
Y,k{LkPLl — j{LlLk,p}), where {.,.} is the anti-commutator 
and the sum is over a set of jump operators L^. We consider 
a dissipator that splits into three parts L = Ll + Ld + Lr. 
Here and Lr describe the coupling to external particle 
reservoirs at the left and right boundaries, respectively, each 
with two jump operators = y/r{l =F f)/2ci^ and LJ^r^ 

^r(l ±/)/2c| ^, where F is the coupling strength, identical 
for both reservoirs, and < / < 1 is the driving bias [29]. We 
consider moderate coupling F = 1 throughout this paper [34]. 
The effect of the driving is to force the system far from equi- 
librium thereby giving a quantum analogue of the well studied 
classical exclusion stochastic model [35, 36]. As depicted in 
Fig. 6 it is composed of two processes; forward pumping cre- 
ating particles at the left edge and removing them at the right, 
and backwards pumping doing the opposite. When / = 
these processes have equal strength and so cancel exactly, re- 
sulting in the stationary solution p = 1 to Eq. (2) with no 
net current, irrespective of A [37]. For / > the bias favours 
forward pumping raising the possibility of a NESS possessing 
a finite current. The remaining contribution Lii accounts for 
bulk dephasing in the chain and is described by a jump opera- 
tor = y/y{l — 2nj) for each site j = 1,...,N with a uniform 
dephasing rate y. 

By directly simulating Eq. (2) and taking the long time 
limit, the state p(f) converges to the time-independent NESS 
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FIG. 2: (a) The current-driving profiles in the weakly-interacting 
regime A = 0.5 for increasing (as indicated by the arrow) dephas- 
ing rates y and N = 16. (b) Identical plot to (a) but in the strongly- 
interacting regime A = 2. (c) The density {n j) of the system when 
A = 2,/= l,andA'= 16, corresponding to y= 0,0.05, 1.00. (d) The 
current as a function of y for driving biases / = 0.1,0.2,0.3 and 0.4 
(from bottom to top curve) with A = 2, and A' = 40. The dashed lines 
correspond to {J)a=o the non-interacting analytic result [44, 45]. 



of the system. A solution can be computed efficiently in a con- 
trolled way and accounting for significant many-body correla- 
tions by applying the time evolving block decimation (TEBD) 
algorithm [38, 39] to a matrix product operator description 
of p(f). This highly compact representation enables the effi- 
cient evaluation of relevant expectation values and makes ac- 
cessible much larger system sizes than exact diagonalization 
or Monte Carlo approaches [40-42]. Moreover TEBD can be 
applied effectively over a large parameter range allowing us to 
examine the properties of the system as a function of / beyond 
the /<C 1 hnear response regime [43]. The transport proper- 
ties are found by computing the current crossing site j from 
the operator Jj — ;(cjcy+i — h.c). In the NESS the current 
(Jj) ~ (J) is homogeneous throughout the system. 

Dephasing enhanced transport — It is known that in the 
weakly interacting regime |A| < 1, in the absence of dephas- 
ing, the system is an ideal ballistic conductor for any driving 
/, with a nearly flat density {nj) and a current (/) oc f which is 
independent of [29, 46]. The introduction of dephasing has 
been shown to induce diffusive transport, where instead 
features a constant gradient and {J) f /N, characteristic of 
an Ohmic conductor [32, 44]. In either case the maximum 
current through the chain occurs at maximal bias f —I, where 
only forward pumping is present, as might be intuitively ex- 
pected. In the non-interacting limit A = it has been proven 
rigorously that a homogeneous chain cannot exhibit any de- 
phasing enhanced transport [44, 47] ; this behaviour was also 
suggested for weakly-interacting systems where |A| < 1 [32]. 
In Fig. 7(a) a complete picture of the latter result is presented. 
Specifically, for A ~ 0.5 the reported current-driving profiles 
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show that dephasing monotonically degrades the current for 
any driving, confirming that this behaviour persists even in the 
presence of weak interactions. Thus, the present work focuses 
on the strongly interacting regime |A| > 1. 

In the absence of dephasing, and for weak driving / <C 1, 
transport was found [46] to be diffusive when |A| > 1, a con- 
troversial finding given that the integrability of the system was 
initially conjectured to lead to ballistic transport [48]. How- 
ever, for strong driving / ^ 1 it was discovered [29, 49] that 
the NESS exhibits a particle domain at the left edge of the 
chain, e.g. dominated by the configuration 1 1 1 • • • 10 • • • 00), 
irrespective of the sign of A. This strongly suppresses the cur- 
rent as (J) oc exp{~N), characteristic of an insulator Conse- 
quently, the current (J) at y = exhibits non-linear behaviour 
with the driving /, leading to an effect known as negative 
differential conductivity (NDC) where increasing the driving 
eventually decreases the current [29, 49]. In Fig. 7(b) the y = 
curve shows that this causes a near complete suppression at 
/ = I for A = 2. Consequently, in the strongly interacting 
regime the system presents the intriguing property that more 
current flows at an intermediate bias / < 1 where some back- 
ward pumping is present. 

The main result of this work is that for |A| > 1 the pres- 
ence of a small bulk dephasing can significantly enhance the 
current. This striking behaviour is illustrated for A = 2 in 
Fig. 7(b) where dephasing up to a moderate rate y w 0.5 is seen 
to increase the current. This elevation in the current is shown 
to occur /or any / > 0, however, it is not uniform in / result- 
ing in the current-driving profile changing with y. Specifically, 
around y « 0.3 the NDC effect is lost and further increases in 
y yield a linear profile in /, resembling a diffusive conductor 
(see Supplementary Material). This indicates that the most 
dramatic enhancement of the current occurs for strong driving 
biases / ^ 1, where the system changes from being insulat- 
ing at y = to yielding the maximal current once y > 0.3. In 
Fig. 7(c) this change at / = I is shown to be coincident with 
the breakdown of the particle domain at the left boundary into 
a near linear density profile even for small dephasing rates. In 
this regime dephasing therefore induces not just quantitative 
increases in the current, but rather causes a major qualitative 
change in the behaviour of the system. 

For a given / there is an optimal dephasing rate yopt, beyond 
which the current is degraded by dephasing. This is illustrated 
in Fig. 7(d), where it is shown that both the maximum en- 
hancement of the current and yopt increase with /. While for 
a weak driving / = 0.1 the current enhancement at its yopt is 
37% relative to y = 0, for / = 1 the enhancement is w 10^. 
For y > yopt the current is reduced because the Ld contribu- 
tion to Eq. (2) dominates over the coherent hopping terms and 
freezes out the dynamics by the Zeno effect [33]. In fact for 
increasingly large dephasing rates the interaction becomes in- 
creasingly irrelevant, as demonstrated in Fig. 7(d) by the con- 
vergence of the NESS current to the exact A = solution with 
dephasing {J)a=o = -2//[(r/4) + (4/r) + {N-l)j] [44, 45]. 

Enhancement mechanism — The origin of the dephasing 
enhanced transport in the strongly-interacting regime is inti- 




FIG. 3: (a) The toy model composed of states | 1) , • • • , | A") coupled 
coherently to their neighbours and with incoherent driving, via an 
auxiliary state |^), between states | 1) o \K). Crucially state \K) 
is elevated in energy by A. See Supplementary Material for precise 
definitions of the model, (b) The current-driving profiles of the toy- 
model as a function of y with A = 2 and K = 20. 



mately connected to the appearance of an insulating NESS at 
/ = 1. This is most easily revealed by examining the very 
strongly-interacting limit |A| ^ 1 with totally biased driving 
/ = 1. Consistent with the numerical findings above, the 
NESS in this extreme limit is insulating for two reasons: (i) 
when |A| 1 there exists an eigenstate | ^D{n)) of H for each 
number sector n, which is a bound state energetically sepa- 
rated from other eigenstates; as « — > A^/2 it becomes exponen- 
tially close, with increasing A^, to being a dark state of / = 1 
driving; (ii) the ergodic nature of the driving and coherent dy- 
namics results in the unique NESS being well approximated 
by a statistical mixture of p = YJ!!=qP,i \ |, with 

the probability p„ exponentially peaked atn=N/2 half-filled 
domain. The formation of an insulating NESS thus arises due 
to an interplay between the nature of the boundary driving and 
the eigenstructure of the chain Hamiltonian. In the Supple- 
mentary Material we show (i) and (ii) by examining the struc- 
ture of the driving and coherent hybridisation of the n parti- 
cle boundary configuration |B„) = |11-10 - - 00) that forms 
|1'd(«)). Crucially, bulk dephasing induces energy dissipa- 
tion at a rate dEy/dt = —2jY,j{c]cj+i +h.c.) proportional to 
the kinetic energy of the state. As a result dephasing incoher- 
ently bridges the spectral gap between \^D{n)) and mobile 
bands of scattering states, thereby dramatically enhancing the 
current as population escapes the dark state. 

The effect outlined is captured concretely by a simple toy- 
model, depicted in Fig. 8(a). The model comprises of states 
|l),|2),...,|/r— l),|/ir), for some size K, with coherent cou- 
pling between neighbouring states. The state \K) is elevated 
in energy by A above the rest, analogous to the configura- 
tion \B„); the other states correspond to breakaway config- 
urations like 1 11 • • • 1010- • -00). This yields a mobile band 
of K — I eigenstates and a single bound state of the form 
I^Vd) w ifjo \2A\-''\K-k), irrespective of the sign of A. 
Via an auxiliary state \s), modelling another particle number 
sector, incoherent driving is included between states 1 1) and 
\K) with an identical / dependence as the full system (see 
Supplementary Material for details). Dephasing is also intro- 
duced which scrambles the phase between | K) and all other 
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FIG. 4: For / = 0. 1 and A' = 40: (a) The half-chain entropy 5 as a 
function of A for different y's, and (b) The optimal dephasing rate 
Yopt as a function of A. The circles indicate TEBD results, and the 

solid line corresponds to the fitted function Yopt °^ (A — Aq)'^, with 
P = 0.819 and A() = 1.07. The inset plots show the generic behaviour 
of (J) with Y above and below Aq. (c) The scaling of Aq with N along 
with the fitted power-law shown where a = 3.906, b = 1.066 and 
c = 0.995. (d) The current for several staggered potential strengths 
B with A = 0.5 (upper panel) and A = 2.0 (lower panel). 



states. Now since the amplitude for 1 1) in \^d) scales as 
|2A| it is exponentially close to being a dark state of / = 1 
driving. Moreover, as / = 1 driving pumps into the majority 
component \K) of | ^o), it leads to an insulating NESS at the 
extremal point / = 1 and y = 0. The current-driving profile 
is shown as a function of / and y in Fig. 8(b), and is seen 
to reproduce both NDC and dephasing enhancement of the 
strongly-interacting regime of the full system. The toy-model 
thus reveals that this behaviour is entirely a consequence of 
the insulator at / = 1 and y = 0. Moreover, it also shows 
that same physics underpins the NDC effect as dephasing en- 
hancement. When / = 1 and y > dephasing induces inco- 
herent transitions to the mobile band. Similarly for y = and 
/ < 1 the leading order effect of the weak backward pumping 
process \s) <— \K), followed rapidly by the stronger forward 
pumping process \s) — > | /T), is to decohere | I'd). 

Signature of a non-equilibrium phase transition — The toy- 
model is only valid for |A| ^ 1 and predicts NDC for any 
|A| > 0. In the non-interacting A = limit the exact solu- 
tion {J)a=o indicates that dephasing acts only to degrade the 
current. However, numerical results presented in Fig. 7(a) 
already indicate that dephasing-degrading behaviour persists 
with weak interactions [32]. This suggests that between the 
regimes |A| ^ 1 and |A| ^ 1 a cross-over in the response of 
the system to dephasing emerges. The point |A| = 1 is of par- 
ticular significance since at zero-temperature it separates the 
equilibrium phases of the underlying spinless fermion Hamil- 
tonian, marking the opening of a gap above the ground state, 
and more generally divides a continuous eigen-spectrum for 



|A| < 1 from one with numerous gaps for |A| > 1 [50, 51]. In 
the driven system it was previously observed that |A| sa 1 also 
separates ballistic and diffusive transport regimes for weak 
driving [29], however directly linking these equilibrium and 
non-equilibrium properties can be misleading [53]. 

For y = 0, studies have attempted to identify the non- 
equilibrium regimes by examining two-point correlations, 
such as Cij = {njnj) — {ni){nj). The strongly-interacting 
regime was found to retain correlations extending over a 
fraction of the entire system as its size increases, while the 
weakly-interacting regime does not, with a sharp cross-over 
found at A^ ~ 0.91 [52]. Here we refine this signature by 
adopting a more general measure of correlations. Specifically, 
we compute the entropy S = — Y,a log2 of '^^e Schmidt 
coefficients Xa arising when the full NESS density operator p 
is factorized into two half-chains as p = La^a<3^<3^ [38, 46, 
55]. The entropy S quantifies both the quantum and classical 
correlations between the two halves of the chain and is read- 
ily accessible within the matrix product operator description 
used in TEBD numerics. For zero dephasing Fig. 9(a) shows 
that S features a peak at A w 1, indicating that a significant 
elevation in correlations occurs as the NESS reorganises itself 
across the expected non-equilibrium phase transition between 
ballistic and diffusive transport there. 

For finite dephasing the entropy S monotonically decreases 
with y, but interestingly, the peak is maintained for moderate 
rates shifting slightly to larger A. Similar to Cjj this distin- 
guishing feature in S is progressively washed out as the sys- 
tem becomes diffusive for all A with dephasing [32, 44, 45] 
(see also Supplementary Material). Despite this, dephasing in 
fact offers an alternative and clear signature of this underly- 
ing non-equilibrium phase transition. In Fig. 9(b) we show 
the optimal dephasing rate yopt for weak driving as a function 
of A, for = 40. A threshold of Aq ~ 1 .07 is apparent where 
for A < Aq, Yopt = 0, indicating dephasing-degrading transport, 
while for A > Aq, Yopt is non-zero and increases monotonically 
with A, indicating dephasing-enhanced transport. The value of 
Ao is size-dependent, and a scaling analysis shown in Fig. 9(c) 
demonstrates that, to good approximation, the isotropic point 
A = 1 separates these regimes in the thermodynamic limit. 
The abrupt on-set of dephasing enhancement therefore occurs 
when the many-body strong correlations are largest. 

Conclusions — We have presented a study of the effects 
of dephasing on the transport properties in a boundary driven 
one-dimensional interacting spinless fermion chain. The ap- 
pearance of a cooperative many-body NESS exhibiting NDC 
provides a previously unexplored form of dephasing enhanced 
transport, distinct in origin from earlier examples in non- 
interacting systems. Using a toy model we isolated the min- 
imum requirements for observing NDC and dephasing en- 
hanced transport in a strongly-interacting system. The ef- 
fect is also unrelated to integrability [56], suggesting that our 
findings it will also apply in related, more realistic models 
of strongly-correlated electrons such as the t — J or Hubbard 
models, and in heat transport [54]. Moreover, tantalising signs 
of dephasing energy dissipation from bound states predicted 
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here have already been observed in a recent cold-atom ex- 
periment [11]. Finally, the prospects of verifying more di- 
rectly both NDC and dephasing enhancement with cold-atoms 
seems promising [26, 27]. 
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Supplementary Material: 
"Dephasing enhanced transport in 
non-equilibrium strongly-correlated quantum 
systems" 

In these appendices we consider in more detail a number 
of points discussed in the main text. Specifically, in Sec. we 
provide evidence that for a strongly interacting regime | A| ^ 1 
at maximal bias / = 1 the transport becomes diffusive in the 
large dephasing limit y > 1. This is followed in Sec. with an 
analysis of the structure of driving at / = 1 between configu- 
ration states. Using this we demonstrate that a dark state exists 
in the limit |A| ^ 1 and we characterise it and the properties of 
the non-equilibrium stationary state (NESS) using perturba- 
tive arguments. In Sec. the essential features of the strongly 
interacting limit are reduced down into a toy-model which is 
seen to exhibit a current-driving profile very similar to the full 
spinless fermion system. The simplicity of the model allows 
an analytic study of its properties and complements the re- 
sults given in the main text. Finally, in Sec. the effects of 
dephasing on the two-point correlations in the full system are 
outlined. 
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FIG. 5: The cuiTent (7) at / = 1, divided by the resulting density dif- 
ference An between the ends of the system (excluding the boundary), 
is plotted against the system size A' — 3 up to A' = 100 sites. The data 
is for dephasing rates Y= 0.5 (x) and Y= 1.0 (■). The solid lines 
show the fitting to (J) /An = k{N — 3)^" for each y. These yield 
K = 1 .288 and a = 0.863 for y = 0.5, and K = 1 .228 and a = 0.958 
fory= 1.0. 



DIFFUSIVE TRANSPORT WHEN |A| > 1 WITH / = 1 AND 

y>0 

In the strongly-interacting regime dephasing is seen to 
modify the current-driving profile from an NDC curve to be- 
ing linear in /. At / = 1 this is also coincident with the break- 
down of the particle domain pinned to the left boundary giving 
way to a linear density profile {rij). Both these features sug- 
gest that dephasing alters the strongly-interacting maximally 
biased regime from being insulating to diffusive. Here we 
present scaling evidence to prove this claim. 

A diffusive current fulfills the usual diffusion equation 
(/) = KV(ny), where K is the conductivity. Given that when 
y > we observe that {tij) is linear in j everywhere but the 
boundary sites, its gradient simplifies to 



V(«,- 



{flN-l) - («2) 



An 

N-3' 



(1) 



where An is the density difference between opposite ends of 
the system (after discarding the boundary sites). This indi- 
cates that a signature of diffusive transport will be a scaling 



(A 

An 



1 



N-3 



(2) 



In Fig. 5 results are presented for the {J) j An at / = 1 with 
A = 2 as a function of the system size up to = 100 sites with 
dephasing rates y = 0.5 and y = 1.0. Also plotted are power- 
law fittings for each y which are seen to accurately model the 
data. These show that the decay of current with A^ is sub- 
diffusive, decaying slower than 1/(A — 3), for y = 0.5 but has 
approached diffusive behaviour once y = 1.0. 



EXISTENCE OF AN APPROXIMATE DARK-STATE AT 
/ = 1 WHEN |A| > 1 

Formally a dark-state | \|/) of an open quantum system is a 
pure state which is simultaneously an eigenstate of the Hamil- 
tonian // 1 \(/) = £ | \|/) and a zero-eigenvalue eigenstate of all 
the jump operators comprising the dissipator i^(|v|/) (v(/|) = 
describing noise acting on the system. Here we show that deep 
in the strongly-interacting limit |A| ::g> 1, with fully polarised 
driving / = 1 as shown in Fig. 6(a), there exists a state | ^d) 
which, as the system size grows, becomes exponentially 
close to satisfying these requirements and is thus an approxi- 
mate dark state. 

In the extreme |A| — > °a limit configuration states such as 
1 10110 - ••Oil), where the particle occupancy on each site of 
the chain is explicitly specified, are exact eigenstates of H. A 
key property of boundary driving is that for any value of /, it 
only incoherently connects configurations within a quadraplet 
of states 1 0x0) , 1 0x1) , 1 1x0) , 1 1x1), where x is any length 
N — 2 occupancy bit string. Thus, if x has [n — l) 1 's then the 
driving couples states within the total number sectors n—l,n, 
and n + 1. This structure constrains the evolution caused by 
the driving processes to shuffling population between states in 
these isolated quadraplets. As depicted in Fig. 6(b), at / = 1 
within each quadruplet, defined by x, there is one configu- 
ration 1 1x0) which, owing to the occupancy on the left and 
empty site on the right, is entirely decoupled from the driv- 
ing, while also being the sink for all driving transitions. As a 
result the effect of this incoherent evolution alone is to even- 
tually drive all the population among the quadraplet of states 
into this dark-configuration. 

Of particular relevance here are the dark-configurations 
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FIG. 6: (a) A schematic of the system, showing the coherent hopping x, the nearest-neighbour density-density interaction A, dephasing y, and 
the incoherent boundary driving in the strong driving limit / = 1 where the bias is maximal. In this case the driving pumps particles into the 
system at site 1 and ejects them at site A', both at a rate F. (b) For any / the driving alone incoherently connects a quadraplet of configurations 
1 0x0) , 1 0x1) , 1 1x0) , I 1x1) . The situation for / = I is shown where the driving can be seen to only pump into the so-called dark configuration 
I 1x0) . (c) For N = 4 sites and / = 1 the complete set of driving quadraplets and the transitions induced by coherent hopping between them from 
the half-filled domain state | 1100) is shown. The vertical axis displays the particle number sectors and energy gaps between configurations. 
We observe that for N = 4 one hop connects | 1100) to | 1010), which is a dark configuration of its quadraplet, and a further hop is needed 
before it is coherently connected to a configuration that is not dark to the / = 1 driving. This disconnection between the half-filled domain 
configuration and the driving only occurs for yV > 4, and so this is the smallest size which displays NDC. 



possessing an n particle domain pinned to the left boundary 
as 

77 

\B„) = I lTr^'nTr QOQ- - -000 ). (3) 

^-77 

For each particle number sector n the configuration | B„) is 
separated from other relevant configurations by an energy gap 
(9(|A|), akin to a domain binding energy. As we move to 
the limit |A| ^ 1, hopping between configuration states oc- 
curs resulting in | B„) giving rise to an eigenstate | ^/)(«)). In 
Fig. 7(a) the exact spectrum for = 12 chain with A = 10 
is shown. For each number sector n the spectrum displays a 
sequence of narrow bands, split by the small finite hopping, 
and separated by gaps of order A. The highest lying eigen- 
state I ^£)(n)) for each sector is also highlighted. Its properties 
are readily determined by treating hopping as a perturbation. 
Specifically, to lowest-order in |2A|^' hopping hybridizes, 
across an energy gap of A, the state \B„) with the configura- 
tion |111---110100- - -000), where the outermost particle of 
the domain has broken away. In Fig. 6(c) the pattern of inco- 
herent driving transitions and coherent hopping for = 4 is 



illustrated. We note that in this case 1 82) = \ 1 100) has no di- 
rect coherent transition to any configurations which couples to 
the / = 1 driving. More generally, so long as A^ > 4 and « > 2 
we observe that the break-away configuration for a domain 
state I B„) is also a dark configuration of its own quadraplet. 

This approach can be taken further by computing higher- 
order corrections due to hopping. Since |A| 1 it is instruc- 
tive to do this approximately by focusing on states of the form 

\nj,k) =cjcl\B„), (4) 

where 1 < j <n and n <k<N, which describe a single break- 
away particle or hole propagating from the domain wall. Since 
the repeated action of hopping on | Z?„) originates around the 
domain wall, and is also detuned by the gap A, we find that 
hopping mixes in configurations where the particle and/or 
hole have hopped a times in total, with an amplitude scaling as 
(9(|2A|^'^). This indicates that the particle/hole propagation is 
suppressed with its distance from the domain wall through the 
unit-filled/empty regions. For each n, the eigenstate 
is therefore predicted by the particle-hole picture to have a do- 
main wall that remains exponentially localized at site «, within 



a length scale ^ ^ l/ln(|2A|) [29], and a deviation 5„(j) from 
the perfect domain configuration | B„) given by 



(a) 



8f(i) 



2(n-j+l) 



2U-") 



, i<y<n 

n<j<N 



(5) 



The validity of this picture is established by comparing this 
to the actual density deviation 8„(j) of the exact eigenstate 
In Fig. 7(b) this is done for \ ^>D{N/2)) with = 
12, and the agreement between 8„(j) and ?>f!^{j) is seen to be 
excellent everywhere but the boundaries. 

Since hopping predominately connects \B„) with dark- 
configurations of the form 1 1x0), the leading order contribu- 
tion to the amplitude in | for a configuration which is 
not dark is (9(|A|-™"(" l"-'^/2|)). This corresponds to whether 
the hole or particle has the shortest path to the left or right 
boundary, respectively. In the |A| 1 limit we therefore 
conclude that the eigenstates with < « < N/2, 
form a hierarchy of states with «, characterised by a decreas- 
ing amplitude for any configuration coupling to the / = 1 
boundary driving. Eigenstates with a domain size 
n scaling with thus become exponentially close to being 
zero eigenstates of the / = 1 driving with increasing system 
size. Of these states the one with n = N/2, where the domain 
spans exactly half the chain, has the most suppressed coupling 
0{\A\^'^^^) to the driving and is thus the closest approxima- 
tion of all of them to an exact dark state. 

Next, we examine the expected structure of the non- 
equilibrium stationary state (NESS) p in the |A| ^ 1 limit. 
The open dynamics of particle number conserving systems 
like that considered here have been studied extensively [29, 
37] with regard to their ergodic and mixing properties estab- 
lishing that a unique NESS exists for any /. The key property 
ensuring this is that for a finite hopping amplitude every con- 
figuration within each particle number sector can be reached 
from any other, while the incoherent ejection/injection of par- 
ticles by the driving connects neighbouring sectors. As a re- 
sult the complete state space of the system can be accessed. 
Furthermore, the NESS for this open system will be block di- 
agonal in the number sectors. At / = 1 the approximate dark 
states I for each sector « are expected to play a promi- 

nent role due to their ability to trap population. This can be 
better understood by approximating the NESS as a statistical 
mixture with probabilities p„ of these eigenstates in each sec- 
tor as 
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FIG. 7: (a) The energy eigen-spectrum of the spinless fermion 
Hamiltonian governing the chain for A' = 12 and A = 10. Ener- 
gies Ev have been shifted by ^A(A'— 1) so that the state | 'i'£i(0)) = 
1 00 . . . 00) has zero energy. The spectrum includes contributions 
from all particle number sectors « = 0, . . . , 12 with each composed of 
sequences of narrow bands, split by the small finite hopping, and sep- 
arated by gaps of order A (since the spectra arising from = 0, • • • ,5 
are identical to those of n = 7, . . . , 12, only the former are shown 
in addition to « = 6). The highest lying eigenstate for each sector 
I (")) is highlighted and seen to be isolated by a gap A from eigen- 
states composed of breakaway configurations \n,i,i). (b) For the 
eigenstate 1 1*0(6)) the exact density deviation 8„(/) from the corre- 
sponding boundary domain configuration \B„) is shown (■), along 
with predicted by considering single particle/hole propagation 

(solid line). The inset shows a schematic of the expected exponential 
localization of the domain wall. 



:X/,„|^o(«))(^o(n) 

«=0 



(6) 



The purity of the NESS in this approximation, tr(p^) = 
I^)i'=o/''n' is reduced only through mixing between sectors. 
The probabilities p„ are then determined by demanding 
that at stationarity there is detailed balance condition be- 
tween the incoherent transition rates connecting neighbour- 
ing number sectors (see inset of Fig. 8). For sector «, as- 



sumed to be frozen in the state |1'o(«)), the output tran- 
sition rates scale with the probability of a hole being at 
the left boundary as F (*d(«) | cic\ \ 1'd(«)) - r|2A|-2", and 
the probability of a particle being at the right boundary as 
r(^D(«) I ^c/v !>!'£,(«)) - r|2A|-2(A'-«). Considering sec- 
tors « — 1 , « and « + 1 we then have the equality of incoming 



9 



and outgoing transitions in n as 

2(N-n) 
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These equations are solved inwards from the extremal n = 
and n = sectors, where \'Vo{Q)) = |00...00) and 
|1'fl(A?)) = |ll...ll),togive 



Pn 
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i2A| 



2\n-N/2Y 



(8) 



for « = 0, 1, . . . and where p = Pn/i is fixed by the nor- 
malization condition Y!!l=QPn = 1- In Fig. 8(a) these pre- 
dicted probabilities of occupation for each sector n are plot- 
ted against the exact values for the NESS with N = 6 and 
A = 10, and found to yield excellent agreement aside from 
the extremal sectors. Owing to the hierarchy of approximate 
dark states this indicates that the NESS will be predominately 
a mixture of | ^o{n)) peaked around the "best" dark state with 
n = N/2. This result also predicts that the purity of p scales 
as 



(9) 



independent of A^. In Fig. 8(b) this prediction is plotted against 
the exact value of 1 — tr(p^) of the NESS for = 6 as a func- 
tion of A, again showing excellent agreement even as A ^ 1 . 

While the discussion of the driven interacting system con- 
sidered has revolved mainly around its NESS properties, the 
existence of approximate dark states at / = 1 will also heavily 
influence its dynamical behaviour far from stationarity. In par- 
ticular, if we envisage starting at time f = in with an empty 
chain, since the growth of a domain requires the propagation 
of holes toward the left boundaries, the increasing suppression 
of this process in each sector « implies that the convergence 
to the NESS will be exponentially slow. This is the physical 
reason why our numerical calculation of the NESS for f — I 
is limited to N ^ 16. Indeed effective suppression of the cur- 
rent sets in once a small domain, « w 5 has formed [29]. As 
a result dephasing enhancement, identical to that seen here in 
the NESS, equally applies to transient dynamical states where 
the domain has not yet saturated at n = N/2. 



THE TOY-MODEL 

The analysis so far has focused on |A| 1 and / = 1. 
In this case we saw that numerous approximate dark states, 
whose occupation is favoured by / = 1 boundary driving, 
cause the stationary state to become insulating. Remaining 
in the strongly interacting limit we now wish to isolate the 
effects on the transport for all drivings /, i.e. the current 
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FIG. 8: (a) The probability p„ on a logarithmic scale of the NESS p 
occupying the n particle sector is shown for an exact calculation (o) 
for A' = 6 and the predicted form (solid line) given in Eq. (8). (b) The 
purity 1 — tr(p-) against A of the NESS p on a logarithmic scale for 
an exact calculation ( x ) for N = 6 and the predicted form (solid line) 
given in Eq. (9). The inset shows a schematic of the detailed balance 
approximation applied to compute these predictions. 



driving profile, due to the existence of such bound states. 
To do so we construct a simple toy model capturing the es- 
sential physics. The structure of the model is motivated by 
considering the half-filled domain state \Bj^/2) and its corre- 
sponding break-away configurations \N/2,i,j) which even- 
tually connect it to the boundary driving. As such the toy 
model is composed of K configuration states \ l) ,\2) , . . . ,\K) 
for some size K >2. To mimic the interaction binding energy 
of |BAr/2), we distinguish the configuration \ K) by elevating 
it in energy by A above the set of otherwise degenerate con- 
figurations \1) ,\2) , . . . ,\K — I) which model the break-away 
states \N/2,i,j). In addition to these configuration energies, 
the state 1 1) , 1 2) , . . . , | A') are also coherently coupled to their 
neighbours via "hopping" processes given by 



K-l 



//, = ! (fc+l|+h.c.). 
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The current operator for the model then follows as 

-/=-'LV) (A:+l|-h.c.), 

k=l 

which measures the flow within the coherently connected con- 
figurations \1) ,\2) ,...,\K). To model the driving in the full 
system, which incoherently connects one particle number sec- 
tor to another, we introduce an auxiliary state | s) whose func- 
tion is simply to be an intermediary. The jump operators de- 
scribing the driving then take the form 

L± = v/r(iT/)/2Aj 4 = v/r(i±/)/2A± 

where = (1 1 and A+ = \K) {s\, with A+ = (A^)''' and 
A^ = (A^)^. Thus, via \s), the driving incoherently induces 
transitions between the boundary configurations 1 1 ) and | K) 
with a bias /. At / = driving in both directions is equal and 
it is easily confirmed that the NESS is p = 1/{K + 1), yield- 
ing (J) = as in the full spinless fermion chain. At the op- 
posite limit, / = 1, population is asymmetrically driven from 
\l) ^ \K). To complete the analogy with the full system the 
toy model also includes dephasing, at a rate y, via the jump 
operator Lz = v^Az where 

Az = l-2\K) {K\, 

whose action is to scramble the phase of any superpositions 
between | K) and the other configurations. A schematic of the 
toy model illustrating all the coherent and incoherent contri- 
butions is shown in Fig. 9(a). 

The current (J) = tr(yp) for the NESS p can be solved ana- 
lytically for the toy model as a function of /, y and A, although 
the complete expression is lengthy. Since the model was mo- 
tivated by the perturbative limit |A| ^ 1 of the full system the 
physically relevant part of this result is found by keeping only 
the lowest order terms in A^ ' . This gives 

(/:-i)(8y/+(i-/)/r) /^y 

In Fig. 9(b) the current-driving profile (/), rescaled by A^, is 
plotted for K ~ 20. Two key features emerge from this result. 
For y = the (1 — /)/ in the numerator of Eq. (10), which 
enforces zero current at the / = and / = 1, causes (J) to 
display negative differential conductivity (NDC). Moreover, 
as the expansion involves only even powers of A^^ this be- 
haviour is independent on the sign of A. Dephasing enhance- 
ment of the (J) is evident from the linear y term in Eq. (10) 
which eventually destroys the NDC effect. However, since 
y and F appear only linearly it indicates that this lowest or- 
der expression is also only valid for weak dephasing and driv- 
ing rates. In particular the dephasing degrading behaviour ex- 
pected for large y, due to the Zeno effect, is not described by 
Eq. (10). The expression is also only valid for K >3 because 
NDC is not seen for K = 2. This is similar to how NDC is only 
seen for > 4 in the full system. Just like in Fig. 6(c) for the 




FIG. 9: (a) The schematic of the toy model and its various processes 
and properties. These include the nearest-neighbour coherent hop- 
ping between the set of states | 1) , 1 2) , . . . , | A"), an energy offset A 
for state \K), incoherent transitions 1 1) o \s) and |A^) -(-> \s) be- 
tween the boundary states and the auxiliary state and dephasing 
at a rate y. The rates for the incoherent driving transitions ^ ( 1 ± /) 
are also listed, (b) The NESS current {J), rescaled by A^, for K = 20, 
as a function of the driving / and dephasing y, described by the ap- 
proximation in Eq. ( 1 0) . The behaviour at / = 1 with y is emphasised 
by the additional (red) line. 

full system, having no direct coherent coupling between the 
boundary configurations where driving occurs, i.e. 1 1) and 
1 3) for = 3, is essential for NDC to emerge. In Fig. 3(b) 
of the main article the exact current-driving profile for the toy 
model is shown for A = 2 as a function of moderate y's, be- 
yond the applicabiUty of Eq. (10). This confirms the wider 
similarity of the response of the toy model to that observed in 
the full spinless-fermion system. 

Given their similar behaviour, the toy model provides a 
tractable means of unravelling the origins of NDC and dephas- 
ing enhancement in the full many-body system. In Fig. 10(a) 
we plot the energy eigen-spectrum of the toy model. By con- 
struction we see that it mimics some of the features seen in 
Fig. 7(a) for a single particle number sector of the spinless 
fermion system. Specifically, there is a high-lying eigenstate 
I I'd), separated by a gap (9(|A|) from a dense band of eigen- 
states. In Fig. 10(b) this band of eigenstates are seen to be de- 
localized over the bulk of the system excluding the boundary 
configuration | K) . For a given size K the eigenstate \ ^d) has 
the form |^d) ~ ifjo \2A\-'' \K - k), to within 0{\2A\--'^), 
as seen in Fig. 10(b). As K ^ °o such an exponential wave 
function is simply the discrete analogue of the well known 
bound state of a ID 8-potential. Given that the amplitude for 
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FIG. 10: (a) The spectrum Ev of the Hamiltonian for the toy model 
for K = 20 and A = 2. Eigenstate v = is the highlighted auxil- 
iary state I s) whose energy has been set arbitrarily. The highest- 
lying eigenstate v = A" is also highlighted and is the approximate 
dark-state of the model | 'Vd)- Beneath this state, separated by a gap 
0(|A|), is a band of eigenstates | ^y) split by the finite hopping, (b) 
The probability distributions | 'Pv ) |^ of the eigenstates over the 
configurations \k) are shown. The band of eigenstates are seen to 
be delocalized over the configurations | 1) , 1 2) , . . . , | A' — 1) and ex- 
punged from the boundary configuration | K) , while characteristic of 
a bound state | "Pfl) is predominately pealced at | A^) with exponential 
tail into bulk. 



the left boundary configuration 1 1) scales as |2A|'^''-, we see 
that I ^o) becomes exponentially close, with increasing K, to 
being a dark state of the driving when / = 1 . 

The driving at / = 1 exclusively pumps into the configura- 
tion I A") whose dominant overlap is with | ^o) . Consequently 
so long as Y = population gets progressively trapped in this 
dark state giving rise to (/) = 0, characteristic of an insulat- 
ing NESS. Remaining at / = 1 and switching on a non-zero 
dephasing directly decoheres the exponentially decaying su- 
perposition within \^d)- This is equivalent to the coherent 
trapping of population, caused by the energetic gap, being by- 
passed by dephasing induced incoherent transitions connect- 
ing I I'o) directly to the delocalized band of eigenstates. Cur- 
rent flow in the system is thus made possible via the ensuing 
non-stationary mixture of these eigenstates. Further increases 
in dephasing eventually degrades the current once the mono- 
tonically decreasing mobility of the delocalized eigenstates, 
caused by the Zeno effect, outweighs the flux of population 
escaping from | ^'o) ■ Since the toy model only has a single 
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FIG. 11: (a) Identical to Fig. 9(b) the NESS current (J), rescaled by 
A^, for = 20 is shown as a function of the driving / and dephasing 
Y zoomed in around the limit / = 1, Y = insulating point. The 
behaviour at / = 1 with Y is emphasised by the additional (red) line, 
(b) In the strong driving limit / ~ 1, the process | A") — > 1 5) occurs 
at a rate oc ( 1 — /) <C 1 making it slow (as indicated by being faded 
out), while the A^ process | ^) — >■ | A") occurs at a rate <x f I and so 
is rapid. The dominant driving process is therefore \K) ^\s) ^\K), 
which has the effect of destroying any coherence | A^) has with other 
states, such as | A" — 1) . 



approximate dark state | ^o) it confirms that its existence, at 
one isolated point / = 1 and Y = 0, is alone enough to make 
the current-driving profile for < / < 1 and Y > exhibit both 
NDC and dephasing enhancement. 

Another key insight from the toy model is that the emer- 
gence of a non-zero current from the insulating point / = 1 
and Y = 0, involves identical physics either when Y is increased 
slightly from zero, or when / is reduced slightly from unity. 
Examining Eq. (10) at / = 1 shows that (J) = 2{K — 
while for Y = an expansion about / = 1 the current is 
(J) ^ j{K — — /)r to lowest order. This suggests a cor- 



respondence 



id 



(11) 



which is confirmed in Fig. 11(a) in an exact calculation for 
K ~ 20. The equivalence is understood by considering the ef- 
fect of the driving slightly below / = 1 where a backward flow 
process is introduced. In this limit the rate of driving from 
|/:) ^ I i) taking the state out of | 'I'd) is given by i(l -/)r 
and therefore slow. In contrast the rate of driving | ^ | K) 
back is 5(1 +/)r, and therefore rapid. Focusing on the dy- 
namics of these two driving processes alone, as in Fig. 1 1(b), 
we suppose that the initial state of the system over the config- 
urations I i) , I A' — 1 ) and I K) has the form 



/ 

p(0)= Pk-i{0) Pk-iAO) 

\ P*K-l,KiO) Pk(P) 

where there is some coherence between \ K — 1) and | K) , but 
no population initially in \s). Evolving this state according 
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FIG. 12: The NESS density-density correlations \C{r)\ across the 
chain as a function of the fraction distance r about the centre for A = 
2 and / = 0. 1. In (a) a moderate dephasing y = 0.05 was included. 
The correlations are long-ranged and extend over a larger fraction 
with increasing A'. This evidenced by the kink, coming from C{r) 
changing sign, moving to larger r. In (b) a strong dephasing y = 1 
was present and the correlations are short-ranged as shown by their 
diminishing size and extent with increasing A'. 

only to these driving processes yields solutions 

p.,(r) = i(l-/)r(l-e-r')p^(0), 

PK{t) = e-V(0) + ^(l+/)r(l-e-^')P^(0), 
PK-i{t) = P/f-i(0), 

We therefore find that for the populations the stationary f — > oo 
limit is approached at a rate F to give (oo) = j ( 1 — /)rp a: (0) 
and Pa:(°o) = 5(1 + /)rp/f (0), while the coherence Pa--i,a'(0 
decays to zero at a rate ^(1 — /)r. Now, by considering only 
the dephasing process we instead find that the populations are 
unchanged while the coherence decays as 

pK-l.K{t)=e-^'^PK-l.K{Q). 

Matching of these decoherence rates fits with Eq. (11). We 
therefore conclude that the emergence of a non-zero current 
when reducing / from unity is, to leading order, caused by 
the resulting decoherence of the dark state |^d), identical to 



the effect of dephasing alone. This behaviour with / around 
/ = 1 , combined with (J) = at / = and the continuity of 
(J) with /, is already enough to imply that NDC behaviour 
will be observed in the current-driving profile. Thus NDC and 
dephasing enhancement are underpinned by the same mecha- 
nism. 

CORRELATIONS AND DEPHASING 



Returning to the full spinless fermion system we now con- 
sider the effect of dephasing on two-point correlations across 
the chain. We find that in addition to modifying the NESS 
current and density profiles, dephasing also leaves a signa- 
ture on the correlations which distinguishes between weakly- 
and strongly-interacting regimes. Specifically we compute 
the NESS density-density correlations Cij = — (n;) {nj), 
where sites ; and j are positioned symmetrically around the 
centre of the chain. This is conveniently represented as C(r) 
where r = |; — y|/A^ is the fractional separation of the points 
for the system size A^, so for example C{r = 1) gives the cor- 
relations between the boundaries. In Fig. 12(a) C(r) is plot- 
ted for different sizes in the strongly interacting regime 
with A = 2, / = 0.1 and a moderate dephasing rate y — 0.05. 
Finite correlations are seen to exist even for a large frac- 
tional separation r, and although they decay with r impor- 
tantly these correlations are retained for a growing fraction 
of the chain as A^ increases. This property has been argued, 
in the y = case, to be evidence that the NESS possess gen- 
uine long-range order [32, 44]. Here our results show that 
this long-range order persists even in the presence of moder- 
ate dephasing, and so correlations similar to that at y = exist 
in the dephasing enhanced NESS. The situation for stronger 
dephasing y = 1, shown in Fig. 12(b), is markedly different. 
Now correlations are smaller and diminish faster with r than 
at weaker dephasing rates. Also the correlations which are 
present cover a smaller fraction r (implying a fixed absolute 
size) as the system size A^ grows. The dephasing has therefore 
degraded the correlations to being short-ranged. Moving to 
the weakly interaction regime with A = 0.5 it has previously 
been found [52] that correlations are short-ranged, similar in 
form to Fig. 12(b), even at y = 0. Indeed a cross-over for C(r) 
from short-ranged to long-ranged correlations at y = was 
found at A > 0.91 [52]. Like the current and density profiles, 
we see that also the correlations in the strongly interacting 
regime with large dephasing rates become increasingly simi- 
lar to those of the weakly interacting regime, washing out this 
cross-over. 



